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The nonlinear spin dynamics in antiferromagnetic crystals is studied for the magnetic structures 
similar to that of hematite. For the case when only two magnetization vectors are non-zero and the 
Hamiltonian has an axial symmetry, a reduction to a Hamiltonian system with one degree of freedom 
is performed, based on the corresponding conservation law. The analysis of the phase portraits of 
this system provides tractable analytical and geometric descriptions of the regimes of nonlinear spin 
dynamics in the crystal. 
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o 
o 

I, INTRODUCTION 

o 
O 

The phenomenological method of describing the properties of crystals with magnetic order is based on attributing 
C*") ' classical magnetization vectors Si (i — l,...n) to the magnetic sublattices of a crystal (n is the number of the 
sublattices), 0]. These vectors are subject to exchange interaction between themselves, anisotropic interaction with 
the crystal lattice, and, optionally, Zeeman interaction with external magnetic fields, [J]. The dynamics of these 
magnetizations can be described via the corresponding Hamiltonian system. Generally, the latter is nonlinear and its 
solving presents substantial difficulties. In the regimes when the spin vectors are close to their equilibrium positions, 
the linear approximation can be used, which allows to obtain the spin wave spectra in a relatively straightforward 
manner. However, the general nonlinear dynamics of sublattice magnetizations is also of considerable interest, e.g. 
for nonlinear regimes of magnetic resonance. Hence there is a need for methods of qualitative and quantitative 
i— i, investigating those dynamics in various situations. 

To this end, the analytical tools devised for the study of mechanical Hamiltonian systems can be employed. In the 
present paper we apply an approach based on the reduction of a Hamiltonian system to another Hamiltonian system 
of lower dimensionality. The key idea is the considering of a subalgebra of the Poisson algebra of dynamic variables. If 
iy— ^ a subalgebra contains the Hamiltonian or, more generally, the Hamiltonian depends on the elements of the subalgebra 
and several first integrals, then the elements of the subalgebra are subject to a new Hamiltonian system with a lower 
number of phase variables. It was Poincare who first performed a reduction of a similar type in the three-body 
problem (see [4|]). A reduction based on considering a Poisson subalgebra was applied to the magnetization dynamics 
' in superfluid 3 He B in (Q), lowering the phase space dimensionality from 6 to 3. 
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o 
o 



In the present work we consider the class of Hamiltonians having an axial symmetry. As a main example we use 
the case described in the classical work by Dzyaloshinsky, - the Hamiltonian comprising the leading terms in 
the magnetic energy of the four sublattices of the antiferromagnetic crystal a-i^C^. It belongs to rhombohedral 
system, having a third-order symmetry axis, The magnetic energy comprises exchange terms and several spin- 
lattice interaction terms, some of them corresponding to the so-called Dzyaloshinsky-Moriya field leading to week 
ferromagnetism of a-Fe^O^ (other possibilities include weak additional antiferromagnetism - the case of CV2O3 in 
and the distortion of magnetic structure described in Q). In our case, if only second-order terms and the largest 
fourth-order term are taken into account (see Q, section "Ferromagnetism of a-FeiOs"), then the Hamiltonian is 
invariant with respect to the rotation about the axis of the crystal. This results into the sum of the longitude angles 
of the spins being a cyclic variable, i.e. not entering the Hamiltonian (for a proper choice of dynamic variables). This 
allows to perform a reduction to a Poisson subalgebra as outlined above, decreasing the number of phase variables by 
2 units. In the case of a-Fe^O^, where only the ferromagnetic vector and one antiferromagnetic vector are non-zero, 
this reduction leads to a system with one degree of freedom, which can be effectively investigated by means of phase 
portraits. That provides a detailed picture of the nonlinear regimes of the spin dynamics in the given approximation. 
The same results hold for the carbonates of Fe, Mn and Co with a difference only in the values of the parameters of 
the Hamiltonian. 
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II. THE HAMILTONIAN FORMULATION 

We follow the well-known paper Q by Dzyaloshinsky for constructing the Hamiltonian system describing the spin 
dynamics in a-Fe^O^. The unit cell of the crystal is rhombohedral. The crystal has 4 magnetic sublattices, and the 
four corresponding Fe ions in the unit cell lie on the body diagonal of the rhombohedron - the third order symmetry 
axis of the crystal. The magnetizations of the sublattices are denoted si,s 2 ,S3,S4. The Poisson brackets between 
their components have the usual form of the brackets for angular momentum: 

{Sia, ~ E-af3^Siy ) {Sia, 0, Z ^ J, (1) 

where the Greek letters denote the Cartesian coordinates of the vectors (here and below repeated indices imply 
summation over 1,2,3). The degeneracy of these brackets leads to the existence of four Casimir functions sf, i — 
1,...,4. 

Following Dzyaloshinsky, we describe the magnetic ordering in terms of the following four vectors: 

Sl + S*2 + S3 + 34 , 
Si - S 2 - S 3 + S*4, 
si - S 2 + S3 - S4, 

si + s 2 - s 3 - s 4 . 

The vector m, the total magnetic moment, corresponds to ferromagnetism (if it is the only non-zero vector, the 
ordering is ferromagnetic), the vectors U - to antiferromagnetism, each of them describing a particular pattern of 
antiferromagnetic ordering, Q. The Poisson brackets for these variables, following from (Q]), read: 

At the temperatures close to that of the antiferromagnetic transition the magnetic energy can be expanded in powers 
of m, h, as their values are small. Let us use the rectangular coordinate system with the z-axis directed along the 
axis of the crystal, the x-axis - along one of the second-order symmetry axes. 

The form of the possible terms in the magnetic energy can be deduced from the analysis of the magnetic symmetry. 
A thorough description of this approach is given in p[. As is shown in p}, in the case of a-Fe^O?, symmetry restrictions 
lead to the following general form of the second-order terms in the magnetic energy of the system: 

E = Ail x + A 2 l 2 + A 3 l% + Brh 2 + ail\ z + a 2 l\ z + a 3 lj z + bm 2 z + (3i(h x m y - h y m x ) + fc{hxh v ~ hyhx)- 

Here the first four terms correspond to the exchange interaction, the other terms - to the relativistic spin-lattice 
interaction and the magnetic dipolar interaction. Experimental results show that in a-Fe^O^ the antiferromagnetic 
structure corresponds to the vector l\, pj], which means that in equilibrium approximately si = — s*2 = — S3 = S4. 
There is also weak ferromagnetism described by the vector rh. So, by physical considerations, the system can be 
restricted to only these two vectors, which leads to the Poisson brackets: 

{m Q , mp} = e Qi a 7 m 7 , {m a , lp} = s a pjl 1: {l a , lp} = e a p 1 m 1 , (2) 

and the Hamiltonian, 0|: 

H = 2 + f + f l * + \ m ' + 13 ( lxl7lv ~ lymx) + 1 ^' (3) 

where I = h, the index 1 is omitted in all positions, and normalizing factors are introduced. It should be noted that 
this Hamiltonian has a forth-order term depending on the main magnetic vector I = l\. This model will be the object 
of our consideration. 



III. SYMMETRY AND REDUCTION 



It is easy to notice that the Hamiltonian J3]) is invariant with respect to the rotation about the z-axis, i.e. the 
anisotropy axis. This leads to m z being a first integral. Indeed, calculation shows that rh z = {m Zl H} = 0. Moreover, 
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the Poisson algebra <[2j) has two Casimir functions: 

(m + f) 2 , {m-lf. 
To take advantage of these facts, let us introduce the following variables: 

m + I t> fa — I 

which are simply g = si + stj,, h = s*2 + S3 in terms of the sublattice magnetizations. Their Poisson brackets have the 
usual form for angular momentum: 

{9a, gp} = Zapyg-f, {h a , hp} = £ Q /3 7 /l 7 , {<? a , hp} = 0. 



The structure of the Poisson algebra given above admits of a reduction that substantially simplifies investigating 
the dynamics of the system. The key point is to cast the Poisson brackets for vector dynamical variables in a form 
that relies on the scalar ones. The idea is essentially due to K. Pohlmeyer, who employed it for studying the 
algebra of currents in field theory. In paper [(| the method was used to study the spin dynamics in the B-phase of 
superfluid 3 He in the regime of turned off magnetic field. 

In our case the reduction is obtained through the following system of variables: 



9 = (V ff 2 - fff cos(/? g , yj g 2 - g 2 z s iny g ,g 2 ), 
h= (y/h 2 - hi cos(p h , ^/h 2 - h 2 z smtp h ,h z ). 

The new variables have the following sense: g - the modulus of g, g z - its z-component, and ip g - its longitude angle 
(in xy-plane). 



Thus, we take into account the symmetry of the system and turn to the scalar quantities, necessary for the 
Pohlmeyer reduction. It is worthwhile to note that a similar transformation was used in papers Q, j§] for studying 
chaotic dynamics in two-spin systems. 



Using the reverse transformation: 



9 = \ hi 



gl+gli fg = arctan 



9y_ 

9x 



we obtain the following Poisson brackets: 



Wg,9z} = 1, {g, 9z} = 0, {g, tpg} = 0, 

{<p h ,h,} = i, {h,h z } = o, {h,cp h } = o, 



all the g-h brackets being zero. 



Thus, we have got two advantages. The first one is the lowering of the phase space dimensionality by two units, as g 
and h are Casimir functions, and their values can be fixed. Therefore, we have a system with two pairs of canonically 
conjugate variables: <p g ,g z and iph,h z . The second one is that we can easily employ the fact that m z — g z + h z is a 
first integral. Indeed, another formulation of this is that the system is invariant with respect to the rotation about 
the z-axis. This means that the Hamiltonian should depend only on the difference of the angles ip g — ipt, not on their 
specific values. In fact, calculation of the Hamiltonian after the above transformations gives: 



9; 



'\A 2 ~ h 2 z cos (ip g - iph) - 2g z h z + g 2 + h* 



2 V9 2 ~ 9zV h2 ~ h l cos (fa - fh) + 2 9zh z + g 2 + h 2 



+ i a ( 9z - hz ) 2 + \b(g z + h z ) 2 - 2/3 y/<F=&y/W=tf z sin [tp g - y h )+ 



He 



-2^/g 2 - gl^h 2 - h 2 cos (ip g - <p h ) - 2g z h z + g 2 + h 2 



Applying the canonical transformation: 



u = — — — — , p 



V2 ' ru V2 
ip g + tp h g z + h 



V2 ' 
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we obtain the Hamiltonian: 

H=\A(- cos (V2u) y/2g 2 - ( Pu + p„) V 2 ^ - (p« = P„) 2 + .9 2 + ^ 2 + P 2 - P 2 ) 
+iS (cos (V2u) V2.9 2 - (p„ + p„) 2 x /2/i 2 - {vu - Pv) 2 + .9 2 + ^ 2 - P 2 + P 2 ) 
+ap2 +bpl-[3 sin \/ 2 5 2 - (p„ + p„) V 2 ^ 2 - (p« = P«) 2 

+ ±C (- cos (>/2u) ^-(p u +p„)V2h 2 -( Pu -p„) 2 + 5 2 + ^ 2 + p£ ~ pi) 2 - 



It is easy to see that v does not enter this function. So, as noted above, p v is a first integral, and as the variables 
u,p u have zero Poisson brackets with it, we conclude that the Hamiltonian [4] defines a separate Hamiltonian system 
in variables u,p u only. In other words, after fixing the integrals of motion the function [4] becomes a member of the 
Poisson subalgebra of dynamic variables generated by u,p u and the first integrals, thus defining Hamiltonian dynamics 
within this subalgebra. 

The study of the dynamics of our system then proceeds as follows. Firstly, we study phase pictures of the reduced 
system for u,p u . Secondly, we consider the "lift" of the obtained two-dimensional dynamics to the six-dimensional 
space of the initial spin variables. Results obtained in this way admit of graphic representation. 

The dynamical regimes of the reduced system can be studied by means of phase portraits in the u — p u phase 
plane. A typical portrait is shown in fig. [TJ It contains several fixed points: centers and saddles, and separatrices, 
each of which either connects a pairs of saddles or forms a loop having the origin and the end in the same saddle. 

A typical trajectory going around a center not far from it is shown in figs. El 03 in the spaces of the vectors I, m 
respectively; a trajectory close to a separatrix - in figs. [fl [H It's easy to see that stationary points in the phase 
portrait correspond to horizontal circles in the spaces of the magnetic vectors I, m, i.e. to circular precession of these 
vectors. Other trajectories manifest more complicated behaviour, as can be seen in figs. [6j [7l 

An especially interesting example is a trajectory going very close to a separatrix. This means that for a long period of 
time it approaches a saddle point, manifesting roughly circular precession of magnetizations, but at a certain moment 
it starts moving away from that saddle performing some complicated dynamics before reaching another almost-steady 
regime of circular precession. 
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FIG. 2: A trajectory surrounding a center (elliptic fixed point) in the phase portrait (fig. [TJ, shown in the space of the 
coordinates of the antiferromagnetic vector I. The curve is winding over a 2D-torus. 
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FIG. 4: A trajectory close to a separatrix connecting two saddles in the phase portrait (fig. [TJ, shown in the space of the 
coordinates of the antiferromagnetic vector I. The vector migrates between two horizontal circles corresponding to the two 
saddles. 
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FIG. 5: A trajectory close to a separatrix connecting two saddles in the phase portrait (fig. [TJ, shown in the space of the 
coordinates of the ferromagnetic vector fa. The vector migrates between two circles corresponding to the two saddles. The 
curve is confined to a horizontal plane m z = const. 




FIG. 6: A trajectory taking S-shaped path between several centers and saddles in the phase portrait (fig. [T]), shown in the 
space of the coordinates of the antiferromagnetic vector I. 




FIG. 7: A trajectory taking S-shaped path between several centers and saddles in the phase portrait (fig. [TJ, shown in the 
space of the coordinates of the ferromagnetic vector m. The curve is confined to a horizontal plane m z = const. 



